Adsorption of rare-gas atoms on Cu(lll) and Pb(lll) surfaces by van der 
Waals-corrected Density Functional Theory 
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The DFT/vdW-WF method, recently developed to include the Van der Waals interactions in 
Density Functional Theory (DFT) using the Maximally Localized Wannier functions, is applied to 
the study of the adsorption of rare-gas atoms (Ne, Ar, Kr, and Xe) on the Cu(lll) and Pb(lll) 
surfaces, at three high-symmetry sites. We evaluate the equilibrium binding energies and distances, 
and the induced work-function changes and dipole moments. We find that, for Ne, Ar, and Kr on 
the Cu(lll) surface the different adsorption configurations are characterized by very similar binding 
energies, while the favored adsorption site for Xe on Cu(lll) is on top of a Cu atom, in agreement 
with previous theoretical calculations and experimental findings, and in common with other close- 
packed metal surfaces. Instead, the favored site is always the hollow one on the Pb(lll) surface, 
which therefore represents an interesting system where the investigation of high-coordination sites 
is possible. Moreover, the Pb(lll) substrate is subject, upon rare-gas adsorption, to a significantly 
smaller change in the work function (and to a correspondingly smaller induced dipole moment) 
than Cu(lll). The role of the chosen reference DFT functional and of different Van der Waals 
corrections, and their dependence on different rare-gas adatoms, are also discussed. 

PACS numbers: 



I. INTRODUCTION 

Understanding adsorption processes on solid surfaces is 
essential to design and optimize countless material appli- 
cations, and to interpret, for instance, scattering exper- 
iments and atomic-force microscopy. In particular, the 
adsorption of rare-gas (RG) atoms on metal surfaces is 
prototypical^ for physisorption processes. Basically, the 
weak binding of physisorbed closed electron-shell atoms, 
such as RG atoms, is due to an equilibrium between at- 
tractive, long-range van der Waals (vdW) interactions 
and short-range Pauli repulsion acting between the elec- 
tronic charge densities of the substrate and the adatomsi^ 

Up to now RG adsorption on many close-packed 
metal surfaces, such as Ag(lll), Al(lll), Cu(lll), 
Pd(lll), Pt(lll),.. have been extensively studied both 
experimentally^"— and theoretically^'"— while, to our 
knowledge, Pb has not, but for the experimental mea- 
surements of Ferralis et ali^ and, very recently, the the- 
oretical investigation of Zhang et al.^ who studied the 
tribological properties of Ne and Kr on the Pb(lll) sur- 
face. The Pb surface is important for practical applica- 
tions: for instance, there is considerable interest in the 
frictional (tribological) properties of gases on Pb at low 
temperatures; in particular, Pb is usedi^"— as a ma- 
terial for the electrodes and as adsorption surfaces in 
nanofriction experiments because it is easy to grow a 
very uniform film already at room temperature and to 
remove the surface contaminants deposited over time on 
the electrodes, thanks to its large diffusion coefficient. 
The Pb(lll) surface also exhibits interesting and unusual 
properties: for instance, one striking finding is the dras- 
tic difference between the sliding friction of Ne and Kr 
mono- or multilayers J^^ii 

In principle, due to the non-directional character of 



the vdW interactions that should be the dominant one in 
physisorption processes, surface sites that maximize the 
coordination of the RG adsorbate atom were expected to 
be the preferred ones, so that it was usually assumed that 
the adsorbate occupies the maximally coordinated hollow 
site. This assumption was also based on the expectation 
that the atom in the hollow site would be closer to the 
surface, thus experiencing a more attractive potential; 
behind this is the notion that the repulsive potential at 
the surface is proportional to the atomic charge density 
and the natural assumption is that the charge density is 
highest at the locations of the atoms, thus making the top 
site energetically unfavored. Calculations where the to- 
tal adatom-substrate interaction is described by the sum 
of empirical binary potentials, which are widely used and 
often give reasonable results for adsorption energies, seem 
to confirm this expectation since the highly coordinated 
hollow sites naturally emerge as the preferred adsorption 
sites for the adatoms. However this picture has been 
questioned by many experimental^"— and theoretical2r22 
recent studies which indicate that the actual scenario is 
more complex: in particular, for Xe and Kr a general 
tendency is found^i^rJ^ for adsorption on metallic sur- 
faces in the low-coordination top sites (this behavior was 
attributed^i^ to the delocalization of charge density that 
increases the Pauli repulsion effect at the hollow sites rel- 
ative to the top site and lifts the potential well upwards 
both in energy and height); for Ar the situation seems 
to be less clear:— for instance, comparison of theoretical 
and experimental results^ would suggest that the hollow 
sites is still favored for Ar on Ag(lll). 

The importance of polarization effects to determine the 
favored adsorption sites was pointed out by Da Silva et 
al.r^ who studied the interaction of RG adatoms with 
the Pd(lll) surface: in fact, for instance, for Xe the 
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polarization is larger in the on-top site, i.e. the larger 
induced dipole moment increases the attractive interac- 
tion between Xe and the metal surface. Therefore, the 
dominant mechanisms appear to be polarization-induced 
attraction and site-dependent Pauli repulsion. The lat- 
ter, being weaker for the on-top site, stabilizes on-top 
adsorptioui^ 

In spite of this recent substantial progress, the under- 
standing of the interaction of RGs with metal surfaces is 
not complete yet.— It is not clear, for instance, whether a 
system exists where high-coordinated site are always pre- 
ferred. Moreover, there have been relatively few studies 
of adsorption geometries for the smaller RGs, although 
these are probably better candidates for the observation 
of high-coordination sites, due to their reduced polariz- 
ability with respect to that of Xe or Kr: in fact, the 
considerable mismatch between the lattice constants of 
the smaller RGs and that of most metal surfaces cause 
most commensurate structures to have multiple atoms 
per unit cell, so that the characterization and interpreta- 
tion of such systems is quite complex. 

Density Functional Theory (DFT) is a well-established 
computational approach to study the structural and elec- 
tronic properties of condensed matter systems from first 
principles, and, in particular, to elucidate complex sur- 
face processes such as adsorptions, catalytic reactions, 
and diffusive motions. Although current density func- 
tionals are able to describe quantitatively condensed mat- 
ter systems at much lower computational cost than other 
first principles methods, they faiU^ to properly describe 
dispersion interactions. Dispersion forces originate from 
correlated charge oscillations in separate fragments of 
matter and the most important component is represented 
by the R^^ vdW interaction)^ originating from corre- 
lated instantaneous dipole fluctuations, which plays a 
fundamental role in adsorption processes of fragments 
weakly interacting with a substrate ("physisorbed" ). 

This is clearly the case for the present systems 
which can be divided into well separated fragments (RG 
atoms and the metal substrate) with negligible electron- 
density overlap. The local or scmilocal character of the 
most commonly employed exchange-correlation function- 
al makes DFT methods unable to correctly predict bind- 
ing energies and equilibrium distances within both the 
local density (LDA) and the generalized gradient (GGA) 
approximations!^ As a consequence, the basic results of- 
ten depend, even at a qualitative level, on the adopted 
DFT functional: for instance, in their ab initio study of 
the interaction of RG adatoms with the Pd(lll) surface. 
Da Silva et alM- found that the on-top site preference is 
obtained by the LDA for all RG adatoms, while the GGA 
functionals (in the FEE and PW91 schemes) yield the on- 
top site preference for Xe, Kr, and He adatoms, but the 
hollow site for Ne and Ar. Typically, in many physisorbed 
systems GGAs give only a shallow and flat adsorption 
well at large atom-substrate separations, while the LDA 
binding energy turns out to be not far from the experi- 
mental adsorption energy; however, since it is well known 



that LDA tends to overestimate the binding in systems 
with inhomogeneous electron density (and to underes- 
timate the equilibrium distances), the reasonable perfor- 
mances of LDA must be considered as accidental. There- 
fore, a theoretical approach beyond the DFT-LDA/GGA 
framework, that is able to properly describe vdW effects 
is required to provide more quantitative resultsJ^ 

In the last few years a variety of practical methods 
have been proposed to make DFT calculations able to 
accurately describe vdW effects (for a recent review, see, 
for instance, ref. [21I ). We have investigated by such 
a method the adsorption of RG atoms on the Cu(lll) 
and Pb(lll) surfaces. Cu(lll) has been chosen because 
of the many experimental and theoretical data available 
(especially for Xe-Cu(lll)), which can be compared with 
ours in such a way to validate the present approach; as 
mentioned above, the less studied Pb(lll) surface could 
be interesting because, given the relatively large Pb lat- 
tice constant (and hence nearest-neighbor surface Pb-Pb 
distance) it represents a good candidate for a system 
where RG atoms are preferably adsorbed on hollow sites 
(the lattice constant of Pb is 4.95 A, compared to 4.09 A 
for Ag, 4.05 A for Al, 3.92 A for Ft, 3.89 A for Pd, and 
3.61 A for Cu). 



II. METHOD 

In this study we include vdW effects within a standard 
DFT approach by using the method proposed in refs. [22I - 

(where further details can be found), hereafter referred 
to as DFT/vdW-WF, by introducing an additional term 
in the exchange-correlation functional as originally pro- 
posed by Andersson et al^ to describe the interactions 
between separate fragments. This contribution, which 
effectively accounts for the dispersion forces both in the 
uniform electron gas and separate atom limits, has the 
form : 
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In the above formulas r„; is the distance between the 
two separate fragments n and /, and Pn(r) is the n-th 
fragment electronic density. The cutoff Tc is introduced 
to remove the divergence of the integral, taking into ac- 
count that, at small momentum values, the interaction is 
highly damped.— 

In our approach all the fragment densities are conve- 
niently rewritten in terms of the Maximally Localized 
Wannier Functions (MLWFs), {w„(r)}, i.e. Pnir) = 
u'^(r). The MLWFs can be obtained from the occupied 
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Kohn-Sham orbitals, generated by a standard DFT cal- 
culation, by means of a unitary transformation which 
minimizes the functional^^ 

^ = ^^n=^{< Wn\r^\Wn > - < W„|r|w„ . 
n n 

(3) 

The unitary transformation conserves the total density, 
which is however partitioned into single localized frag- 
ments, each of them being characterized by its spread Sn 
and center of mass position r„. It is therefore possible 
to express the vdW correction (see eqs. ^ and ((U) as 
a sum of single contributions coming from each pair of 
Wannier functions belonging to different fragments, by 
approximating the shape of the n-th Wannier functional 
with an H-like exponential. 

The DFT/vdW-WF method has been already suc- 
cessfully applied to several systems, including small 
molecules, bulk, and surfaces j^^—i^i"— in particular it 
allowed us to study the interaction of Ar with graphite 
and of Ar, He, and H2 with Al surfaces (2^*21 of water with 
the CI- and H-terminated Si(lll) surfaces;^ and of RG 
atoms and water with graphite and graphenei^ 

We here apply the DFT/vdW-WF method to the case 
of adsorption of Ne, Ar, Kr, and Xc atoms on the Cu(lll) 
and Pb(lll) surfaces. All calculations have been per- 
formed with the Quantum-ESPRESSO^S ab initio pack- 
age (MLWFs have been generated as a post-processing 
calculation using the WanT package^). Similarly to 
DaSilva et al.^ we modeled the clean and RG-covered 
metal surfaces using a periodically-repeated hexagonal 
supercell, with a (\/3 x •\/3)i?30° structure and a sur- 
face slab made of 15 metal (Cu or Pb) atoms distributed 
over 5 layers (repeated slabs were separated along the di- 
rection orthogonal to the surface by a vacuum region of 
about 24 A). The Brillouin Zone has been sampled using 
a 6 X 6 X 1 fc-point mesh. In this model system the RG 
coverage is 1/3, i.e. one RG adatom for each 3 metal 
atoms in the topmost surface layer. The (a/3 x -\/3)i?30° 
structure has been indeed observed at low temperature 
by LEED for the case of Xe adsorption on Cu(lll) and 
Pd(lll}^ (actually, this is the simplest commensurate 
structure for RG monolayers on close-packed metal sur- 
faces and the only one for which good experimental data 
exist), and it was adopted in most of the previous ab 
initio studies^"— li^ii^ii^. Since the lateral interactions 
between RG adatoms do not play a critical role in the 
RG adsorption site preference;^ for the sake of simplic- 
ity, we have used the same structure also for the other 
RGs (Ne, Ar, and Kr) and in the case of adsorption on 
Pb(lll) as well. 

The Pb or Cu surface atoms were kept frozen (of course 
after a preliminary relaxation of the outermost layers of 
the clean metal surfaces) and only the vertical coordi- 
nates of the RG atoms, perpendicular to the surface, 
were optimized, this procedure being justified by the fact 
that only minor surface atom displacements are observed 
upon physisorptionj^ii^i^ Moreover, the RG atoms were 
adsorbed on both sides of the slab: in this way the sur- 



face dipole generated by adsorption on the upper surface 
of the slab is cancelled by the dipole appearing on the 
lower surface, thus greatly reducing the spurious dipole- 
dipole interactions between the periodically repeated im- 
ages (previous DFT-based calculations have shown that 
these choices are appropriateiSili) . Note that, appar- 
ently, in their recent study of Ne and Kr on Pb(lll), 
Zhang et alr^ have instead considered adsorption on a 
single side of the metal slab; the effect of such a choice is 
non- negligible: in fact, for instance, in the case of Xe on 
Pb(lll), we find that the (absolute value of the) binding 
energy is reduced by 7 meV (about 4 %) with respect to 
that obtained when Xe is adsorbed on both sides of the 
slab. The results of ref. [l^ are thus likely affected by the 
artificial dipole-dipole interactions discussed above. 

We have carried out calculations for various separa- 
tions of the RG atoms adsorbed on high-symmetry sites, 
namely hollow (on the center of the triangle formed by 
the 3 surface metal atoms contained in the supercell), top 
(on the top of a metal atom), and bridge (intermediate 
between two nearest-neighbor metal atoms). Actually, 
two kinds of hollow sites are present: HCP hollows, char- 
acterized by having atoms directly beneath them in the 
next layer of atoms, and FCC hollows where this con- 
dition does not apply; however the HCP-hollow and the 
FCC-hollow sites can be considered equivalent for ad- 
sorption because of the small differences in the adsorption 
properties (for instance, Righi and Ferrario fusing LDA, 
found a difference of less than 1 meV in the adsorption 
energy and of 0.01 A in the equilibrium distance for RGs 
adsorbed on Cu(lll)). For a better accuracy, as done in 
previous applications on adsorption processes^i^i^^i^ 
we have also included the interactions of the MLWFs 
of the physisorbed fragments not only with the MLWFs 
of the underlying surface, within the reference supercell, 
but also with a sufficient number of periodically-repeated 
surface MLWFs (in any case, given the decay of the 
vdW interactions, the convergence with the number of 
repeated images is rapidly achieved). Electron- ion inter- 
actions were described using norm-conserving pseudopo- 
tentials: in the case of Pb and Cu we have explicitly in- 
cluded 14 and 11 valence electrons per atom, respectively 
(those coming from the Sd^", 6s^, 6p^ atomic orbitals for 
Pb, and 3d^°, 4s^ for Cu). As a reference DFT functional 
we chose PW91— because it is widely used in ab initio 
DFT calculations of solids and surfaces and, in particular, 
was adopted in some previous simulations^^ of Xe inter- 
acting with the Cu(lll) surface, which facilitates com- 
parison with the results of the present calculations (note 
that typically PW91 gives similar results to that obtained 
by PBE,— which represents another popular GGA func- 
tional). Using the PW91 functional in test calculations 
with bulk Pb and Cu, for the equilibrium properties the 
agreement with experimental estimates is comparable to 
that found in previous DFT calculationsj ^'^^i"'^^ 

By generating the MLWFs for the Cu(lll) and 
Pb(lll) substrates, we observe a clear quantitative sep- 
aration between the spreads of the MLWFs describing 
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d-like orbitals and those of the (much more deloeahzed) 
MLWFs describing the s- and p-hke orbitals; moreover, 
given the high valence-electron density, screening effects 
are certainly relevant in these metal surfaces. Therefore, 
at variance with previous calculationS j^^i^^'^^i^^ we have 
applied the DFT/vdW-WF correction by explicitly con- 
sidering only the more localized MLWFs corresponding 
to the d-like orbitals, while the s- and p-like electrons arc 
supposed to give a screening-effect^^ contribution, which 
is taken into account by multiplying the vdW correction 
(the Cg coefficients) by a simple Thomas-Fermi factor: 
Jtf = e^^^^^^'^^^^TF where r^^ is the Thomas-Fermi 
screening length relative to the electronic density of a 
uniform electron gas ("jellium model") equal to the av- 
erage density of the s- and p-like electrons of the present 
systems, Zs is the average vertical position of the topmost 
Cu or Pb atoms, and z is the vertical position, measured 
with respect to Zs, of the adatom. In practice it turns 
out that only the topmost metal layer gives a relevant 
contribution, while the effects of the other ones is dra- 
matically reduced by the exponential factor, in line with 
the common expectation about screening effects in metal 
surfacesi^ This observation can be exploited to consider- 
ably reduce the computational cost of the vdW correction 
since only the topmost MLWFs must be really taken into 
account. 



III. RESULTS AND DISCUSSION 

In Tables I- VI results are reported for all the systems 
under consideration, for adsorption on hollow, top, and 
bridge sites. The binding energy. Eh, is defined as 

Eb = l/2{Etot - [E, + 2Erg)) (4) 

where Es^rg represent the energies of the isolated frag- 
ments (the substrate and the RG atoms) and Etot is the 
energy of the interacting system, including the vdW- 
correction term (the factors 2 and 1/2 are due to the 
adsorption of RG atoms on both sides of the slab). 

One should point out that the experimentally mea- 
sured adsorption energy, Ea, includes not only the in- 
teraction of RG atoms with the substrate but also lat- 
eral, vdW, RG-RG interactions however in most of 
previous calculations the mostly attractive lateral in- 
teraction contribution was not considered. As pointed 
out, for instance, by Lee et al.^ who studied n-butane 
on transition-metal surfaces (another typical weak ph- 
ysisorption system where vdW interaction is the only at- 
tractive force between the nonpolar molecule and the sub- 
strate) lateral adatom-adatom interaction energies can 
be as large as 25% of the total adsorption energy at full 
coverage. Ea is here defined as: 

Ea~Eh + {El — Ef), 

where Ei is the total energy (per atom) of the 2D RG 
lattice (that is as in the adsorption configurations but 
without the substrate) and Ef is the energy of an iso- 
lated RG atom. Clearly the quantity in parenthesis in 



the above formula represents the lateral adatom-adatom 
interaction energy (per atom). Note that, in their DFT 
study of Ne and Kr on Pb(lll), Zhang et alw^ seem in- 
stead to identify Ea with E^. 

Ef, has been evaluated for several adsorbate-substrate 
distances; then the equilibrium distances and the cor- 
responding binding energies have been obtained by fit- 
ting the calculated points with the function: Ae~^^ — 
Cs/iz — zq)^, a, B, C3, and zq being adjustable param- 
eters (as illustrated for the Xe-Cu(lll) and Xe-Pb(lll) 
cases in Figs. 1 and 2). Typical uncertainties in the fit 
are of the order of 0.05 A for the distances and a few 
meV for the minimum binding energies. Our results are 
compared to available theoretical and experimental esti- 
mates and to corresponding data obtained using a pure 
PW91 functional, the simple LDA functional, and the 
"seamless" vdW-DF method of Langreth et alS^ (note 
that the vdW-DF method has been also used in the re- 
cent DFT study of Zhang et alM). As can be seen in 
Figs. 1 and 2, and in Tables I and II, the effect of the 
vdW correction computed by DFT/vdW-WF is a much 
stronger bonding than with a pure PW91 scheme, with 
the formation of a clear minimum in the binding energy 
curve at a shorter equilibrium distance. In spite of the 
clear shortcomings of the pure PW91 scheme, in general 
the preferred adsorption site seems to be correctly deter- 
mined by the latter, although the differences between the 
binding energies of the different adsorption sites are very 
small. 

We have also computed Ea (assuming a full monolayer 
coverage of RGs) in the case of Xe on Cu(lll), where a 
RG over layer in the (VS x V3)R30° structure is experi- 
mentally found^ and in the case of Xe on Pb(lll), where 
the formation of a commensurate Xe monolayer was also 
observed]^ As can be seen in Table III, all the meth- 
ods, but pure PW91, correctly predict a smaller Ea on 
Pb(lll) than on Cu(lll), although the quantitative re- 
sults considerably depend on the adopted scheme: in fact, 
pure PW91 clearly underestimates Ea, DFT/vdW-WF 
and vdW-DF give comparable results, while LDA is close 
to DFT/vdW-WF and vdW-DF for Xe on Cu(lll) but 
underestimates for Xe on Pb(lll): this can be explained 
by the fact that LDA is not able to describe properly the 
lateral interactions of Xe adatoms which are further from 
each other on Pb(lll) than on Gu(lll). 

Concerning the adsorption on the Cu(lll) surface (see 
Table I), all the methods used predict that the top config- 
uration is energetically favored in the case of Xe, while 
for Ne, Ar, and Kr the differences among the binding 
energies of the different adsorption configurations are 
quite small (using vdW-DF the same is true also for Xe); 
since these differences are probably comparable to the 
expected accuracy of the calculations, a precise assign- 
ment of the favored adsorption site is not possible. In 
contrast, the hollow configuration is instead clearly fa- 
vored by all the methods (see Table II) in the case of the 
adsorption on Pb(lll) of all the considered RG atoms 
(actually, with DFT/vdW-WF, for Ar on Pb(lll) the 
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bridge site appears to be lower in energy: however, given 
the small difference with the energy of the hollow site, 
this result should not be overemphasized). Our results 
for Pb(lll) are in qualitative agreement with those of 
Zhang et alJ^ who predict that Ne and Kr indeed prefer 
high-coordination hollow sites. Note that the energy dif- 
ference between the hollow and top sites increases by sub- 
sequently considering the PW91, vdW-DF, DFT/vdW- 
WF, and LDA methods (see also Da Silva et al^). 

Interestingly, in the case where several experimental 
reference values are available, namely Xe on Cu(lll), 
our DFT/vdW-WF method performs better (considering 
both the binding and adsorption energy, and the equilib- 
rium distance, see Tables I, III, and IV) than all the other 
schemes: in fact LDA gives reasonable binding ener- 
gies but underestimates the equilibrium distances, while 
vdW-DF underestimates the binding energies and overes- 
timates the equilibrium distances, in line with the behav- 
ior reported for systems including a metallic surface.— 
Also note that, at a variance with the experimental find- 
ings, the vdW-DF method predicts that the top site (see 
Table I) is only marginally favored (and the distance only 
marginally different) than the hollow ones; in general, for 
all the RG atoms on Cu(lll) vdW-DF gives almost iden- 
tical binding energies for the top and hollow adsorption 
sites. In the case of RGs on Pb(lll) the hollow structure 
is favored also by vdW-DF, although the difference in 
the binding energy with respect to the top site is smaller 
than with the present DFT/vdW-WF scheme (the dif- 
ference was instead larger, see the last column of Table 
II, in the study of Zhang et ali^, who used vdW-DF but 
with a reference DFT functional differing from ours by 
the exchange term). In the case of Ar on Cu(lll) and on 
Pb(lll), we observe that our computed binding energies 
compare favorably with the estimates obtained, using a 
simple Lennard- Jones potential, by Cheng et al.,— who 
predicted a binding energy between 70 and 85 meV for 
Ar on noble metals. 

As expected, wc find that, both for adsorption on 
Cu(lll) and Pb(lll), the binding energy increases by 
going from Ne to Xe, in line with the increasing polariz- 
ability of this atom sequence. In particular, for several 
close-packed transition-metal surfaces the binding energy 
of Xe is foundiS to be about 2 to 3 times larger than that 
of Kr, and Ar, respectively, a behavior which is well re- 
produced by our DF/vdW-WF method (the factors arc 
1.5 and 3, and 1.6 and 2.5, for adsorption on Cu(lll) 
and on Pb(lll), respectively). This general behavior is 
also in line with the results of Zhang et ali^ 

Our energetic results are not far from the "best esti- 
mate" reported by Vidali et al^ for Xe on Cu(lll), i.e. 
a binding energy of -183 ±10 meV at a distance of 3.60 
±0.08A (these values represent averages over different 
theoretical/experimental estimates). In their tables Vi- 
dali et al^ also report for Ar on Cu(lll) a binding energy 
of -85 meV at a distance of 3.53 A and for Kr on Cu(lll) 
a binding energy of -119 meV, in fair agreement with our 
results. Lazic et alt^ studied the adsorption of Xe on 



Cu(lll) by a DFT approach where vdW corrections were 
included using the method of Andersson et al.^ using 
PW91 and PBE as reference DFT functional (see the 
last column in Tables I and IV). As can be seen, our re- 
sults arc much closer to the experimental estimate than 
those of Lazic et al.^ which tend to overestimate the 
binding energy and underestimate the equilibrium dis- 
tance. The Xe-adsorbed Cu(lll) surface has been also 
recently investigated by Sun and Yamauchii^ using DFT 
with semiempirical vdW corrections: they found reason- 
able equilibrium distances, however the computed bind- 
ing energy was very overestimated (it was even larger 
than that obtained by LDA) and the favored adsorption 
site was incorrectly predicted to be the hollow site, prob- 
ably due to the use of semiempirical pair potentials which 
favor close-packed structures and high coordinated sites 
(see discussion above). 

From Tables I, III, and IV, on can also see that the 
binding energies are reasonably reproduced by the LDA 
scheme for RGs on Cu(lll), a behavior common to sev- 
eral physisorption systems. However, as already outlined 
above, this agreement should be considered accidental: 
the well-known LDA overbinding, due to the overesti- 
mate of the long-range part of the exchange contribu- 
tion, somehow mimics the missing vdW interactions; the 
equilibrium distances predicted by LDA are clearly un- 
derestimated since LDA cannot reproduce the be- 
havior in the interaction potential. For RGs on Pb(lll), 
the LDA binding energies are instead underestimated as 
a consequence (as discussed above) of the larger equilib- 
rium distances than for RGs on Cu(lll). 

As already found elsewhere for all the used 
schemes, the binding energies correlate with the RG- 
metal distance: in fact, for a given RG, the configura- 
tions having the strongest binding are characterized by 
the shortest RG-substrate distance. Moreover, all the 
methods predict that Ar and Xe adatoms get closer to 
the Cu(lll) surface when adsorbed on top site, as found 
in some previous studies i^iiSiH Remarkably, this behav- 
ior cannot be reproduced^ii^ using a hard-sphere model, 
indicating that there is a significant interaction between 
the Ar and Xe atoms and the Cu(lll) surface so that a 
simple stacking (hard-sphere) model of weakly or nonin- 
tcracting spheres is not valid (for comparison, in Tables 
IV and V wc also list the sums of the RG atoms and 
metal atom vdW literature radii). Instead, for adsorp- 
tion on Pb(lll), the adatoms in the hollow site arc closer 
to the surface than in the top one, in line with the usual 
behavior. These results can be easily elucidated by ana- 
lyzing the parameters of the adopted fitting function (see 
above), Ae~^^ — C3/{z — zq)^ : we find that, as a general 
rule, at the equilibrium distance, the repulsive potential 
term is weaker on the favored adsorption site (for in- 
stance the top site for Xe on Cu(lll) and the hollow one 
for Xe on Pb(lll)), in agreement with the results of Da 
Silva et al^ 

Ferralis et alJ^ studied the structural and thermal 
properties of Xe on the Pb(lll) surface by LEED. They 
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observed the formation of a Xe monolayer with an incom- 
mensurate hexagonal structure with a lattice parameter 
similar to that found in bulk Xe (4.33 A); this structure 
is aligned with the substrate lattice but has a larger unit 
cell, similarly to the case of Xe on Ag(lll), which is also 
an aligned incommensurate monolayer. They also found 
that the heat of adsorption for the first Xe layer is -191 
±10 meV with an, ovcrlayer-substrate spacing of 3.95 
±0.10 A. Looking at Table III we found that our com- 
puted Ea (-205.5 meV for the hollow adsorption site) 
is close to the experimental value and in better agree- 
ment than with the other methods, although our model 
structure is not exactly the same observed experimen- 
tally; moreover, also the Xe-Pb(lll) distance (3.93 A) 
is in excellent agreement (see Table V) with that esti- 
mated by Ferralis et al.r^ which gives further support 
to the reliability of our DF/vdW-WF method. As ex- 
pected, it has been founcU^ that a hard-sphere model is 
unable to give a good description of adsorption of Xe on 
Pb(lll). For Xe-Pb(lll) the heat of adsorption is lower 
than for Xe on any surface measured so fai^ with the 
possible exception of Al(llO) and for alkali metals; a low 
heat of adsorption is not particularly surprising since the 
Pb atoms are much larger than most other metals (the 
vdW radius of Pb is 2.02 A, compared to 1.72 A for 
Ag, 1.72 A for Pt, 1.63 A for Pd, and 1.40 A for Cu), 
implying that the repulsive Xe-Pb interaction prevents 
the Xe from approaching the deeper part of the attrac- 
tive holding potential. It must be noted that Ferralis et 
alJ^ were unable to determine the preferred adsorption 
site, the lack of satellite intensities in the LEED patterns 
indicating that the overlayer is quite uniform and the 
corrugation is small. 

An important quantity which often provides reveal- 
ing details of the bonding mechanism in adsorption pro- 
cesses is represented by the electron density difference, 
An(r) = njiG/s{^) — "s(r) — n^Q{v), obtained from the 
electron density (at the equilibrium geometry) of the RG 
on the substrate, of the clean substrate, and the isolated 
RG monolayer, respectively. Our approach in this re- 
spect is not fully self-consistent because we use the elec- 
tron density obtained at a pure PW91 level, that is with- 
out vdW corrections, however, the effects due to the lack 
of self-consistency are expected to be negligible because 
the rather weak and diffuse vdW interactions should not 
substantially change the electronic charge distribution 4S 
Plots of An(r) for Xe on Cu(lll) and Xe on Pb(lll), 
both in the hollow and top site (see Figs. 3 and 4), show 
that, in agreement with what found prcviousl}*i2i for RGs 
on Pd(lll), the electron density redistribution is stronger 
on the Cu atoms for the Xe on the top site than for the 
hollow; both sites exhibit a significant depletion of elec- 
tron density centered about the Xe atom together with 
a slight density accumulation close to the center of the 
Xe atom, this effect being attributecti^ to orthogonaliza- 
tion of Xe states to the states of the substrate atoms. 
Moreover, for Xe in the on-top site, there is a significant 
electron density accumulation between the Xe atom and 



the topmost surface layer. Interestingly, there is a clear 
tendency of Xe to induce a much larger charge delocal- 
ization on the Cu(lll) surface than on Pb(lll), in line 
with the delocalization mechanism invoked^ii^ to explain 
the preference for the top adsorption site on Cu(lll). 

Since polarization effects are assumed to play a key 
role in determining the favored adsorption siteSfi^ii^ we 
have also computed the change of the work function, 
AW, of the Cu(lll) and Pb(lll) substrate upon adsorp- 
tion of RG atoms. The work functions have been calcu- 
lated as the difference between the averaged electrostatic 
Coulomb potential at the midpoint of the vacuum region 
of the slab and the Fermi energy for the clean Cu(lll) 
and Pb(lll) surfaces we estimate a work function of 4.85 
and 3.86 eV, respectively, in excellent agreement with the 
reference values, that are in the range from 4.90 to 5.01 
cVi^ for Cu(lll), and 3.83 eV for Pb(lll).i^ AW can be 
related to the dipole moment induced in the substrate by 
the presence of the RG adatom, A/x, using the Helmholtz 
equation)^ 

where ^(ixi) is the area of the (1 x 1) surface unti 
cell (in A^) and 8 is the RG coverage; if AW is given 
in eV, then A/i is in debyes. In our case & = 1/3, so 
that A^ = s/Sal/lGnAW, where gq is the Cu or Pb lat- 
tice constant. Our computed AW and A/i values are 
listed in Table VI. In agreement with previous ab initio 
calculations r'liS we find that the RG adsorption induces 
a decrease in the work function, thus indicating that the 
RG atoms behave as adsorbates with an effective positive 
charge; note that this is consistent with the depletion of 
the electron density about the Xe atom discussed above, 
which corresponds to an induced surface dipole moment 
that points out of the surface. For Xe on Cu(lll) our 
estimated AW and A/i values (see Table VI) agree well 
with the experimental estimates^ of -0.60 eV and -0.24 
D, respectively. As can be seen in Table VI, the absolute 
value of A/i increases from Ne to Xe, because the corre- 
sponding electronic polarizabilities increase, and is larger 
for the optimal adsorption site, for instance the top for Xe 
on Cu(lll) and the hollow for Xe on Pb(lll). Moreover, 
it is considerably larger on Cu(lll) than on Pb(lll) in 
line with the energetic analysis reported above, that in- 
dicated a stronger interaction of RGs with the Cu(lll) 
surface than with Pb(lll). 

Zhang et alr^ explain the much larger mobility of Ne 
overlaycrs on Pb(lll), as observed in friction experi- 
ments, than of Kr ovcrlayers on the basis of the different 
activation energies which characterize the lateral motion 
of Ne and Kr atoms on the Pb(lll) surface. The activa- 
tion energies for a monolayer can be directly calculated 
from the difference in the binding energy of the adatom 
between the favored (hollow) site and the transition state, 
which is expected to correspond to the bridge site Con- 
sidering the differences between the binding energy of 
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FIG. 1: Binding energy of Xe on Cu(lll) in the top and hol- 
low configuration using pure PW9f (full and empty circles, 
respectively) and DFT/vdW-WF (solid and dashed line, re- 
spectively), as a function of the distance z from the surface. 



FIG. 2: Binding energy of Xe on Pb(lll) in the top and hol- 
low configuration using pure PW91 (full and empty circles, 
respectively) and DFT/vdW-WF (solid and dashed line, re- 
spectively), as a function of the distance z from the surface. 



the hollow and bridge configurations for Nc and Kr on 
Pb(lll), we qualitatively confirm the trend observed by 
Zhang et al.r^ being our estimated activation energies 
(1.3 meV for Ne and 6.0 meV for Kr ) of the same order 
of magnitude as those reported in ref. |T6| (0.7 meV for Ne 
and 2.5 mcV for Kr). However, such small energy values 
are comparable to (or even smaller than) the expected 
accuracy of the computed binding energies, thus mak- 
ing quantitative estimates of the hopping probabilities^^ 
(which depend exponentially on the aforementioned ac- 
tivation energies) rather questionable. 



IV. CONCLUSIONS 

In summary, by analyzing the results of our study of 
the adsorption of RG atoms on the Cu(lll) and Pb(lll) 
surfaces, one can conclude that the inclusion of the vdW 
corrections by the DFT/vdW-WF method systematically 
improves upon the estimates for the binding energies as 
obtained by a standard GGA approach. In particular, 
using a pure PW91 functional the binding is underesti- 
mated in all cases, while equilibrium distances are over- 
estimated. For all the system considered the vdW cor- 
rection term represents the dominant part of the bind- 
ing energy, although, particularly for RG adsorption on 
Pb(lll), the pure PW91 approach gives a substantial 
contribution. However, vdW interactions appear not to 
play a critical role in the adsorption site preference (the 
same result has been obtained by Zhang et alJ^ study- 
ing the interaction of Ne and Kr on Pb(lll)): Xe on 
Cu(lll) clearly prefers the top site, while for Ne, Ar, an 
Kr on Cu(lll) the differences in binding energies rela- 
tive to different adsorption sites arc so small that is not 




FIG. 3: Electron density difference of Xe on Cu(lll) in (a) 
hollow and (b) top site shown in a plane perpendicular to the 
surface, within the range of ±1 X 10"*e/bohr^ Red (light 
grey) and blue (dark grey) represent electron accumulation 
and depletion, respectively. The green and orange spheres 
indicate the Xe and Cu atoms, respectively. 
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(a) 




easy to attribute a definitive preference; instead, tlic hol- 
low configuration tends to be preferred for adsorption 
of all the considered RGs on Pb(lll), in agreement with 
previous calculations and experimental observations! ^^'^^ 
Moreover, the Pb(lll) substrate is subject, upon rare- 
gas adsorption, to a significantly smaller change in the 
work function, and to a correspondingly smaller (in abso- 
lute value) induced dipolc moment, than Cu(lll). Given 
these relevant peculiarities of the Pb(lll) surface, where 
the hollow site is undoubtedly favored for adsorption of 
RG atoms, this surface would represent an ideal sub- 
strate to study, both theoretically and experimentally, 
high-coordination adsorption sites. 



(b) 




FIG. 4: Electron density difference of Xe on Pb(lll) in (a) V. ACKNOWLEDGEMENTS 

hollow and (b) top site shown in a plane perpendicular to the 
surface, within the range of ±1 X 10""*e/bohr^ Red (light 
grey) and blue (dark grey) represent electron accumulation 

and depletion, respectively The green and grey spheres indi- ^e thank very much F. Costanzo and F. Toigo for 

cate the Xe and Pb atoms, respectively. useful discussions 
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TABLE IL Binding energy, in meV, of RC atoms on the 
Pb(lll) surface computed using the standard DFT-PW91 
calculation, and including the vdW corrections using our 
DFT/vdW-WF method, compared to the LDA result, the 
vdW-DF method by Langreth et alM- and available theoreti- 
cal and experimental (in parenthesis) reference data. 



system 


PW91 


DFT/vdW-WF 


LDA 


vdW-DF 


ref. 


Ne-Pb(lll) 


hollow 


-31.2 


-59.8 


-49.4 


-71.4 


-51.6" 


Ne-Pb(lll) 


top 


-27.8 


-49.1 


-42.9 


-63.3 


-46.8" 


Ne-Pb(lll) 


bridge 


-19.8 


-58.5 


-49.1 


-64.6 




Ar-Pb(lll) 


hollow 


-23.5 


-82.4 


-78.3 


-100.8 




Ar-Pb(lll) 


top 


-22.1 


-75.0 


-64.2 


-95.3 




Ar-Pb(lll) 


bridge 


-22.7 


-84.5 


-76.6 


-100.1 




Kr-Pb(lll) 


hollow 


-30.8 


-132.8 


-98.8 


-136.9 


-134.9" 


Kr-Pb(lll) 


top 


-29.1 


-109.8 


-81.6 


-130.9 


-125.1" 


Kr-Pb(lll) 


bridge 


-24.0 


-126.8 


-96.7 


-136.1 




Xe-Pb(lll) 


hollow 


-59.6 


-193.5 


-142.0 


-192.2 


-172.6" 


Xe-Pb(lll) 


top 


-56.3 


-186.4 


-116.1 


-186.4 




Xe-Pb(lll) 


bridge 


-52.7 


-188.9 


-138.6 


-191.2 





"reflU 
Refill 



TABLE I: Binding energy, in meV, of RC atoms on the 
Cu(lll) surface computed using the standard DFT-PW91 
calculation, and including the vdW corrections using our 
DFT/vdW-WF method, compared to the LDA result, the 
vdW-DF method by Langreth et al^ and available theoreti- 
cal and experimental (in parenthesis) reference data. 



system 


PW91 


DFT/vdW-WF 


LDA 


vdW-DF 


ref. 


Ne-Cu(lll) 


hollow 


-17.6 


-31.6 


-55.7 


-56.1 






Ne-Cu(lll) 


top 


-17.5 


-31.1 


-55.4 


-55.9 






Ne-Cu(lll) 


bridge 


-17.6 


-31.0 


-55.3 


-56.1 






Ar-Cu(lll) 


hollow 


-13.0 


-67.8 


-88.9 


-106.6 






Ar-Cu(lll) 


top 


-13.0 


-71.9 


-94.5 


-106.3 




-85" 


Ar-Cu(lll) 


bridge 


-13.0 


-70.6 


-89.4 


-106.4 






Kr-Cu(lll) 


hollow 


-20.3 


-134.2 


-117.6 


-135.7 






Kr-Cu(lll) 


top 


-20.3 


-131.1 


-126.0 


-135.8 




-119" 


Kr-Cu(lll) 


bridge 


-20.3 


-130.0 


-118.4 


-135.7 






Xe-Cu(lll) 


hollow 


-22.9 


-194.5 


-199.3 


-167.4 




-276^ -268'= 


Xe-Cu(lll) 


top 


-23.1 


-208.1 


-221.9 


-167.7 


-280^ 


-183", -277^= (-190^=) 


Xe-Cu(lll) 


bridge 


-17.1 


-191.2 


-201.0 


-167.4 




-278*" 



"refH 
''ref p 
'^ref 
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TABLE III: Adsorption energy {Ea, see text for the defini- 
tion), in meV, of Xe atoms on the Cu(lll) and Pb(lll) sur- 
faces computed using the standard DFT-PW91 calculation, 
and including the vdW corrections using our DFT/vdW-WF 
method, compared to the LDA result, the vdW-DF method 
by Langreth et al^ and available experimental (in parenthe- 
sis) reference data. 



system 


PW91 


DFT/vdW-WF 


LDA 


vdW-DF 


ref. 


Xe-Cu(lll) hollow 
Xe-Cu(lll) top 
Xe-Cu(lll) bridge 


-51.4 
-51.6 
-45.6 


-289.3 
-302.9 
-286.0 


-297.3 
-319.9 
-299.0 


-268.9 
-269.2 
-268.9 


(-227") 


Xe-Pb(lll) hollow 
Xe-Pb(lll) top 
Xe-Pb(lll) bridge 


-62.5 
-59.2 
-55.6 


-205.5 
-198.4 
-200.9 


-147.9 
-122.0 
-146.9 


-252.2 
-246.4 
-251.2 


(-191") 



"refill 



TABLE IV: Equilibrium RG adatom-surface distance, in A, 
on the Cu(lll) surface computed using the standard DFT- 
PW91 calculation, and including the vdW corrections using 
our DFT/vdW-WF method, compared to the LDA result, the 
vdW-DF method by Langreth et al^ and available theoreti- 
cal and experimental (in parenthesis) reference data; the sum, 
s, of the vdW radii of the RG atom and the Cu atom is also 
reported. 



system 


PW91 


DFT/vdW-WF 


LDA 


vdW-DF 


ref. 


s 


Ne-Cu(lll) 


hollow 


3.90 


3.59 


3.10 


3.70 




2.94 


Ne-Cu(lll) 


top 


3.90 


3.57 


3.09 


3.68 




2.94 


Ne-Cu(lll) 


bridge 


3.90 


3.60 


3.10 


3.68 




2.94 


Ar-Cu(lll) 


hollow 


4.50 


3.48 


3.19 


3.90 




3.28 


Ar-Cu(lll) 


top 


4.50 


3.45 


3.15 


3.86 


3.53" 


3.28 


Ar-Cu(lll) 


bridge 


4.50 


3.43 


3.19 


3.86 




3.28 


Kr-Cu(lll) 


hollow 


4.50 


3.32 


3.21 


3.99 




3.42 


Kr-Cu(lll) 


top 


4.50 


3.36 


3.17 


3.99 




3.42 


Kr-Cu(lll) 


bridge 


4.50 


3.35 


3.20 


3.99 




3.42 


Xe-Cu(lll) 


hollow 


4.70 


3.42 


3.00 


4.10 


3.40^ 3.3T= 


3.56 


Xe-Cu(lll) 


top 


4.40 


3.36 


2.90 


4.09 


3.45^ 3.2^*, 3.25'= (3.60"=) 


3.56 


Xe-Cu(lll) 


bridge 


4.70 


3.41 


3.00 


4.10 




3.56 



"refll 

''ref 
'^ref 
"^ref 
'^ref 
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TABLE V: Equilibrium RG adatom-surface distance, in A, 
on the Pb(lll) surface computed using the standard DFT- 
PW91 calculation, and including the vdW corrections using 
our DFT/vdW-WF method, compared to the LDA result, the 
vdW-DF method by Langreth et al^ compared to the LDA 
result, and available theoretical and and available theoretical 
and experimental (in parenthesis) reference data; the sum, s, 
of the vdW radii of the RG atom and the Pb atom is also 
reported. 



system 


PW91 


DFT/vdW-WF 


LDA 


vdW-DF 


ref. 


s 


Ne-Pb(lll) hollow 


3.80 


3.41 


3.10 


3.70 


3.5" 


3.56 


Ne-Pb(lll) top 


4.00 


3.68 


3.40 


3.90 


3.8" 


3.56 


Ne-Pb(lll) hndge 


3.80 


3.36 


3.27 


3.50 




3.56 


Ar-Pb(lll) hollow 


4.40 


3.68 


3.40 


4.00 




3.90 


Ar-Pb(lll) top 


4.40 


4.04 


3.60 


4.22 




3.90 


Ar-Pb(lll) bridge 


4.50 


3.77 


3.43 


4.10 




3.90 


Kr-Pb(lll) hollow 


4.40 


3.69 


3.40 


4.14 


3.8" 


4.04 


Kr-Pb(lll) top 


4.40 


3.98 


3.70 


4.24 


3.9" 


4.04 


Kr-Pb(lll) hndge 


4.30 


3.79 


3.51 


4.13 




4.04 


Xe-Pb(lll) hollow 


4.30 


3.93 


3.50 


4.30 


(3.95") 


4.18 


Xe-Pb(lll) top 


4.50 


4.02 


3.70 


4.30 




4.18 


Xe-Pb(lll) hndge 


4.70 


3.93 


3.55 


4.31 




4.18 



"reflU 



TABLE VL Work-function change, in eV, and induced dipole 
moment (in parenthesis), in debyes, for RGs adatoms on the 
Cu(lll) and Pb(lll) surfaces, at equilibrium geometries. 



system 


hollow 


top 


Ne-Cu(lll) 
Ar-Cu(lll) 
Kr-Cu(lll) 
Xe-Cu(lll) 


-0.04 (-0.02) 
-0.28 (-0.13) 
-0.54 (-0.24) 
-0.53 (-0.24) 


-0.03 (-0.01) 
-0.37 (-0.17) 
-0.37 (-0.17) 
-0.57 (-0.26) 


Ne-Pb(lll) 
Ar-Pb(lll) 
Kr-Pb(lll) 
Xe-Pb(lll) 


-0.03 (-0.03) 
-0.10 (-0.08) 
-0.11 (-0.09) 
-0.13 (-0.11) 


-0.03 (-0.03) 
-0.03 (-0.03) 
-0.05 (-0.04) 
-0.04 (-0.03) 



